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We study the inclusive hadrodroduction of D° , D + , D* + , and D$ 



qq 1 to-leading order in the parton model of quantum chromodynamics endowed with 



universal non-perturbative fragmentation functions (FFs) fitted to e + e — annihila- 

\J^) , ti° n data from CERN LEP1. Working in the general-mass variable-flavor-number 

^— ^ , scheme, we resum the large logarithms through the evolution of the FFs and, at the 

-x^, 1 same time, retain the full dependence on the charm-quark mass without additional 

' theoretical assumptions. In this way, the cross section distributions in transverse 

£L| momentum recently measured by the CDF Collaboration in run II at the Fermilab 

' , Tevatron are described within errors. 

Ph. 

1. Introduction 

- 1 — i ■ 

. Recently, there has been much interest in the study of charmed-hadron 

| (X c ) production at hadron colliders, both experimentally and theoretically. 

The CDF Collaboration measured the differential cross sections da/dpx for 
the inclusive production of D°, D + , D* + , and Df mesons (and their an- 
tiparticles) in pp collisions at the Fermilab Tevatron (run II) as functions 
of transverse momentum (jpt) in the central rapidity (y) region. 1 Until 
recently, the most advanced theoretical predictions, 2,3 based on quantum 
chromodynamics (QCD) at next-to-leading order (NLO), consistently un- 
dershot all the D°, D + , and D* + data by significant amounts, as is evident 
from Fig. 3 of Ref. 1, while no predictions for mesons existed. Espe- 
cially in view of future physics at the CERN Large Hadron Collider, where 
the continuum production of charmed hadrons will provide important back- 
grounds for numerous new-physics signals, it is an urgent task to deepen 
our understanding of the inclusive hadroproduction of charmed hadrons on 
the basis of QCD in order to render the theoretical predictions as reliable as 
possible, so as to establish a sturdy anchor for new-physics searches. Here, 
we report on recent progress in this direction. 4 
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2. General-Mass Variable-Flavor-Number Scheme 

We wish to advocate the general-mass variable-flavor-number (GM-VFN) 
scheme, which has recently been elaborated for the photoproduction 5 and 
hadroproduction 6 of heavy-flavored hadrons. In this approach, one starts 
from the pt 3> m c region and absorbs the large logarithms ln(p^/m^) 
into the parton density function (PDF) of the c-quark in the incoming 
hadrons and the fragmentation functions (FFs) for the c — ► X c transitions. 
After factorizing the In to;; terms, the cross section is infrared safe in the 
limit m c — > 0, and n/ = 4 is taken in the strong-coupling constant a s 
and the DGLAP evolution equations. The remaining to c dependence, i.e. 
the m^/ptp power terms, is retained in the hard-scattering cross sections. 
These terms are important in the intermediate pr region, where PT^m c , 
and are expected to improve the precision of the theoretical predictions. 
The large logarithms are absorbed into the PDFs and FFs by subtraction 
of the collincarly (mass) singular terms. However, in order to define a 
unique factorization prescription, one also has to specify non-singular terms. 
This is done by requiring that, in the limit px — *• oo, the known hard- 
scattering cross sections of the zero-mass variable-flavor-number (ZM-VFN) 
scheme are recovered. To achieve this, subtraction terms are derived by 
comparing the fixed-flavor-number (FFN) theory in the limit m c — > with 
the ZM-VFN theory, implemented in the MS factorization scheme. 5 ' 6 This 
matching procedure is useful, since all commonly used c-quark PDFs and 
FFs are defined in the ZM-VFN scheme. The latter can then be used 
consistently together with hard-scattering cross sections calculated in the 
GM-VFN scheme. 

3. Numerical Results 

We are now in a position to present our numerical results for the cross 
sections of inclusive D°, D + , D* + , and D+ hadroproduction to be directly 
compared with the CDF data, 1 which come as distributions da/dpT at cm. 
energy >/S = 1.96 TeV with y integrated over the range \y\ < 1. For each 
X c species, the particle and antiparticle contributions are averaged. We 
work in the GM-VFN scheme with n/ = 4, thus excluding X c hadrons 
from Vfc-hadron decays, which are vetoed in the CDF analysis. 1 We set 
to c = 1.5 GcV and evaluate a^ lf \fj,n), where [ir is the renormalization 
scale, with = 328 MeV, 8 corresponding to a {5 \m z ) = 0.1181. We 
employ proton PDF set CTEQ6.1M from the CTEQ Collaboration 8 and 
the NLO FFs 7 that were recently fitted to LEP1 data taking the starting 
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scales for the DGLAP evolution to be fio = m c , m^. We distinguish between 
the initial- and final-state factorization scales, hf and fj,' F , so that we have 
three unphysical mass scales altogether. Our default choice is ^ — \ip — 
Hp = rriT, where hit = \Jvt + m c i s the transverse mass. In order to 
conservatively estimate the theoretical error due to the scale uncertainty, 
we independently vary the values of fiR/mr, lip /rriT, and /ip/rriT between 
1/2 and 2, and determine the maximum upward and downward deviations 
from our default predictions. 

Our theoretical predictions are compared with the CDF data in Fig. 1. 
The four frames refer to D°, D + , D* + , and Df mesons. In all cases, we 
find good agreement in the sense that the theoretical and experimental er- 
rors overlap, i.e. the notorious discrepancy between experiment and theory 
mentioned in Sec. 1 has disappeared. In fact, our theoretical predictions 
provide the best description of the CDF data obtained so far. 

4. Conclusions 

In conclusion, the GM-VFN scheme resums large logarithms by the DGLAP 
evolution of non-perturbative FFs, guarantees the universality of the latter 
as in the ZM-VFN scheme, and simultaneously retains the m c -dependent 
terms of the FFN scheme without additional theoretical assumptions. 
Adopting this framework in combination with new fits of D°, D + , 
and D+ FFs 7 to LEP1 data, we managed for the first time to reconcile the 
CDF data on the production of these mesons in Tevatron run II 1 with QCD 
within errors and thus eliminated a worrisome discrepancy. Furthermore, 
we presented the first NLO predictions for the Df data. 1 
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Figure 1. Comparison of the CDF data 1 
D°,D+,D*+,Df. The solid lines represent < 
lines indicate the unphysical-scale uncertainty. 
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with our NLO predictions for X c = 
■ur default predictions, while the dashed 
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